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(57) Abstract: Some embodiments of the invention provide a location -determination system tliat includes several transmitters and 
at least one receiver. Each transmitter transmits a signal that includes a unique periodically -repeating component, and the receiver 
receives a reference signal. Based on ihe received reference signal, the lucaliun-determinatiun system identifies an estimated location 
of the receiver as follows. For each transmitter in a set of transmitters, the system computes a phase offset between the received 
reference signal and a replica of the transmitter's periodically-repeating component. The system also identifies an a^jroximate 
location of the receiver and an approximate receive time for the received .signal. 'ITic .sy.stcm then uses the identified approximate 
location and lime, and the computed phase oftsets, lu compute pseudonmges for the set of transmitters. Finally, the system identifies 
the estimated location of the receiver by using the computed pseudoranges. 
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Location-Determination Method and Apparatus 

FIELP OF THE mVENTION 

The invention is related to location-determination method and apparatus. 

BACKGROUND OF THE INVENTION 

In many applications, it is necessary to estimate the location of objects in their 
enviromnent. To date, numerous location-determination systems have been proposed 
for this task. One such system is the global positioning system (GPS). This system 
includes a number of satellites that orbit the Earth. It also includes GPS receivers, 
monitoring stations, and differential GPS receivers on Earth. 

GPS satellites transmit signals jfrom which GPS receivers can estimate their 
locations on Earth. A GPS satellite signal typically includes a composition of: 
*(1) carrier signals, (2) pseudorandom noise (PRN) codes, and (3) navigation data. 
(jPS satellites transmit at two carrier frequencies. The first carrier frequency is 
approximately 1575.42 MHz, while the second is approximately 1227.60 MHz. The 
second carrier frequency is predominantly used for military applications. 

Each satellite uses two PRN codes to modulate the first carrier signal. The 
first code is a coarse acquisition (C/A) code. This code is a sequence of 1023 
elements, called chips^ and it modulates at an approximately IMHz rate. The second 
code is a precise (P) code, which repeats on a seven-day cycle and modulates at a 
lOMHz rate. Different PRN codes are assigned to different satellites in order to 
distinguish GPS signals transmitted by different satellites. 
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The navigation data is superimposed on the first carrier signal and the PRN 
codes. The navigation data is transmitted as a sequence of five firames. Each fiame is 
six seconds long, and it contains a time stamp for when the frame was transmitted. 
The navigation-data fiames also provide information about the sateUite's clock errors, 
the satellite's orbit (i.e,, ephemeris data) and other system status data, A GPS satellite 
receives its ephemeris data firom monitoring stations that monitor ephemeris errors in 
its altitude, position, and speed. 

Based on the signals transmitted by the GPS satellites, GPS techniques 
typically estimate the location of a GPS receiver by using a triangulation method, 
which typically requires the acquisition and tracking of at least four satellite signals at 
the 1.57542 GHz fi-equency. GPS signal acquisition often involves correlation 
calculations between the received GPS signal and the C/A code of each satellite at 
various phase offsets and Doppler-shifted firequencies. Code phase and phase delay 
are other names for phase offset. For each satellite, GPS acqmsition techniques record 
the largest-calculated correlation value as well as the phase offset and Doppler-shifted 
frequency resulting in this value. GPS acquisition techniques then identify several 
satellites (e.g., the four satellites) that resulted in the highest-recorded correlation 
values. 

After signal acquisition, a signal tracking process decodes the sigoals fix>m the 
identified sateUites at the phase offsets and Doppler values associated with the 
recorded correlation values for these satellites. Specifically, the signal tracking 
process uses the identified phase offset and Doppler values for the identified satellites 
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to extract each identified satellite's navigation data. Part of the extracted data is the 
time stamp information. 

From the extracted time data, the signal tracking process can compute the 
distance between the receiver and the identified sateUites. In particular, a satellite's 
signal-transmission delay (/.e., the time for a signal to travel jftom the satellite to the 
receiver) can be calculated by subtracting the satellite's transmission time the 
satellite's extracted time stamp) from the time the receiver received the satellite's 
signal. In turn, the distance between the receiver and a satellite can be computed by 
multiplying the satellite's signal-transmission delay by the speed of light TTie 
estimated or exact distance between the receiver and a satellite is often referred to as 
the satelUte's pseudorange. 

After signal tracking, a triangulation process typically computes the location of 
the GPS receiver based on the computed pseudoranges and the locations of the 
satellites. The location of each satellite identified during signal acquisition can be 
calculated from the satellite's ephemeris data. Theoretically, triangulation requires the 
computation of pseudoranges and locations of only three sateUites. However, 
triangulation methods often use the pseudoranges and locations of four sateUites 
because of inaccuracies in the receiver clock. 

Some GPS systems also improve flieir accuracy by using diflferential GPS 
techniques. Such a technique requires the operation of differential GPS receivers at 
known locations. Unlike regular GPS receivers that use timing signals to calculate 
their positions, the differential GPS receivers use thefr known locations to calculate 
timing errors due to the signal path. These differential GPS receivers determine what 
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the travel time of the GPS signals should be, and compare them with what they 
actually are. Based on these comparisons, the differential GPS receivers generate 
"error correction" factors, which they relay to nearby GPS receivers. The GPS 
receivers then factor these errors into their calculation of the transmission delay. 

5 Signal tracking has a number of disadvantages. For instance, it is 

computationally intensive, and hence quite time consuming. Also, signal tracking at 
worst needs 6 seconds, and on average needs 3 seconds, of GPS data to extract a 
selected satellite's time stamp, as the time stamp is embedded in the navigation-data 
frames that are 6 seconds long. This requirement, in turn, reduces the speed of the 
10 location-determination process. Finally, signal tracking and reliable decoding of data 
require higher signal power than that needed for acquisition. This becomes a 
significant impediment when the GPS signal is attenuated, as in indoor or urban 
environments. 

Therefore, there is a need for a global positioning system that can quickly 
15 identify a GPS receiver's location. There is also a need for a global positioning system 
that only needs a small amount of received GPS data to identify the GPS receiver's 
location. More generally, there is a need for a location-determination system that 
addresses some or all of the above-mentioned needs. 
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SUMMARY OF THE INVENTION 

Some embodiments of the invention provide a location-determination system 
that includes several transmitters and at least one receiver. Each transmitter transmits 
a signal that includes a unique periodically-repeating component, and the receiver 
5 receives a reference signal. Based on the received reference signal, the location- 
determination system identifies an estimated location of the receiver as follows. For 
each transmitter in a set of transmitters, the system computes a phase offset between 
the received reference signal and a replica of the transmitter's periodically-repeating 
component. The system also identifies an approximate location of the receiver and an 
1 0 approximate receive time for the received signal. The system then uses the identified 
approximate location and time, and the computed phase offsets, to compute 
pseudoranges for the set of transmitters. Finally, the system identifies the estimated 
location of the receiver by using the computed pseudoranges. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The novel featuies of the mvention are set forth m the appended claims. 
However, for purpose of explanation, several embodiments of the invention are set 
forth in the following figures. 

5 Figure 1 illustrates a signal processing circuit that receives a GPS signal and 

generates a number of digital samples Scorn this GPS signal. 

Figure 2 illustrates a process that determines the location of the GPS receiver 
fi*om the digital samples generated by the signal-processing circuit. 

Figure 3 illustrates a location-determination system. 

1 0 Figure 4 illustrates a process that calculates the approximate location of GPS 

satellites at a particular transmission time. 

Figure 5 illustrates one satelUte that is "over" an approximate location and 
another satellite that is not "over" the approximate location. 

Figure 6 illustrates a process that identifies satelUtes that are over an 
1 5 approximate location. 

Figure 7 illustrates a process that computes phase offsets for several satellites. 

Figure 8 illustrates the ^e of code segments that a GPS receiver receives 
firom several GPS satellites. 

Figure 9 illustrates a process for computmg pseudoranges. 
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Figure 10 illustrates a PRN code received by a GPS receiver, a replica of tbis 
code stored by the GPS receiver, and a phase-shifted version of this replica. 
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DETAILED DESCRIPTION OF THE INVENTION 

The invention provides location-detennination method and apparatus. In the 
following description, numeroxis details are set forth for purpose of explanation. 
However, one of ordinary skill in the art will realize that the invention may be 
5 practiced without the use of these specific details. For instance, some embodimoits of 
the invention are described below by reference to global positioning systems. One of 
ordinary skill will understand tiiat other embodiments of the invention are used in 
other types of location-determination systems. In other instances, well-known 
structures and devices are shown in block diagram form in order not to obscure the 
10 description of the invention with unnecessary detail. 

Some embodiments of the invention are location determination systems that 
determine location of receivers from "reference signals" received by the receivers. As 
used in this document, a reference signal means any type of signal from which 
location information may be derived. Hence, the reference signal can be a GPS 
1 5 ("global positioning system") signal, a CDMA ("code division multiple access") 

signal, a GSM ("global system for mobile communication") signal, a TDMA ("time 
division multiple access") signal, a CDPD ("cellular digital packet data") signal, or 
any other signal fcom Miiich location information may be derived 

In the embodiments described below, however, the reference signal is a GPS 
20 signal that can be used to estimate die location of GPS receivers. On Earth, a GPS 
receiver typically receives a GPS signal that is a composite of several signals 
transmitted by GPS satellites that orbit the Earth. The characteristics of such GPS- 
satellite signals were described in the background section. 

8 
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The GPS receiver can be a standalone device, can be part of another mobile 
device (e.g., a personal digital assistant ("PDA"), wireless telephone, etc.), or can 
communicatively connect to another mobile device (e.g., connect to a Handspring 
Visor PDA through its proprietary Springboard). Several such architectures for the 
5 GPS receiver are described in United States Patent Application, entitled "Method and 
Apparatus for Determining Location Using a Thin-Client Device," filed on December 
4, 2000, and having serial number 09/730,324. The disclosure of this application 
United Stated Patent Application, entitled "Method and Apparatus for Determining 
Location Using a Thin-Client Device," filed on December 4, 2000, and having serial 
10 number 09/730,324) is incorporated herein by reference. 

Some embodiments estimate the location of a GPS receiver by (1) initially 
generating K-digital samples Sj, sk firom the GPS reference signal received by the 
receiver, and then (2) using the digitized GPS reference data to estimate the location 
of the GPS receiver. The GPS receiver typically performs the digitization operation. 
15 Also, in some embodiments, the GPS receiver digitizes only a portion of the received 
GPS signal. An example of a signal-processing circuit that a GPS receiver can use to 
generate such samples will be described below by reference to Figure 1. 

Some embodiments estimate the location of the GPS receive firom the 
digitized GPS reference data by using a location-determination process that can be 
20 performed either (1 ) completely by the GPS receiver, (2) completely by another device 
or computer in communication with the GPS receiver, or (3) partially by tiie GPS 
receiver and partially by another device or computer in communication with the GPS 
receiver. 
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In some mibodiments, the locadon-detennination process computes 
pseudoranges for a set of satellites without performing signal tracking on the received 
GPS reference data. In some of these embodiments, the process computes for each 
satellite in the set a phase offset between the samples of the received reference signal 
5 and a replica of the satellite's PRN code. The system also identifies an approximate 
location of the receiver and an approximate receive time (e.g.^ an approximate time 
for the start of the received signal). 

The system then uses the receiver's approximate location, approximate receive 
time, and the computed phase offsets to compute pseudoranges for the set of satellites, 

10 To compute a satellite's pseudorange, some embodiments use the satellite's computed 
phase offset and the receiver's approximate location to calculate the length of all or 
most of the PRN codes transmitted by the satellite from the time that it transmitted the 
signal resulting in the first sample 5/ to the time the receiver generated the first sample 
sj. Finally, the system identifies the estimated location of the receiver by using the 

1 5 computed pseudoranges. 

In certain circumstances, the estimated receiver location matches the exact 
receiver location. In other circumstances, the estimated receiver location matches the 
exact receiver location to such a high degree of accuracy that it is indistinguishable 
from the exact location to an observer. In yet other situations, however, the estimated 
20 location differs from the actual location of the GPS receiver by a certain error amount; 
in these situations, some embodiments take steps to ensure that this error (between the 
estimated and actual receiver locations) is tolerable for the particular location- 
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determination application. Several more specific location-determination processes 
will be explained by reference to Figures 2-10. 

I. GENEEIATING THE DIGITAL REFERENCE DATA 

Figure 1 illustrates a signal processing circuit 100 that receives a GPS signal 
and generates K-digital samples from the GPS signal. As shown in this figure, the 
signal processing circuit 100 mclxides a GPS antenna 105, a GPS tuner 1 10, a clock 
1 15, a down-converter 120, and an analog to digital ("A/D") converter 125. 

The GPS antenna 105 receives a GPS signal s(t), which on Earth is a 
composite of noise and several signals transmitted by several GPS satellites that orbit 
the Earth. The antenna 105 and its associated circuitry are configured to receive tiie 
reference GPS signal s( t) at a GPS carrier frequency, which currently is around 1 .57 
gigahertz ("Ghz"). The RF tuner 110 receives the GPS signal s(t) from the antenna 
105. This tuner 1 10 is tuned to capture signals at the approximate frequency of the 
GPS signal. Hence, the tuner captures the GPS reference signal s(t) received by the 
antenna 105. 

The RF tuner communicatively couples to the clock 1 15 to receive a clock 
signal. The clock 115 generates one or more clock signals to synchronize the 
operation of the components of the signal-processing circuit 100. This clock also 
receives a synchronizmg clock signal 130. The clock 1 15 can receive the 
synchronizmg signal 1 30 from a variety of sources. For instance, in some 
embodiments, the signal-processing circuit 100 includes an RF processing circuit tiiat 
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(1) captures a radio signal with the synchronixdng signal, and (2) supplies this signal to 
the clock 115. 

The synchroni2dng signal allows the clock to set initially its internal time. It 
also allows the clock to synchronize with a reference clock. In some embodiments, 
5 the receiver is part of a network (e.g. , a wireless telephone network) and its clock 115 
is synchronized with the network clock so as to keep the receiver clock within a few 
micro- or milli-seconds off the GPS clock. Some embodiments synchronize the clock 
1 15 to be within 1 millisecond or less of the GPS clock. The degree of the receiver 
clock's inaccuracy will depend on (1) source and manner of obtaining the 
10 synchronizing signal 130, (2) the accuracy of the clock of the synchronizing-signal 
source, and (3) the accuracy of the GPS-receiver clock. 

It should be noted that other embodiments use other techniques to maintain the 
receiver clock within a desired time interval of GPS time. For instance, some 
embodiments might have the receiver estimate and compensate for its clock drift 
15 internally. Other embodiments use high-precision clocks in their receivers. 

As shown in Figure 1, the down converter 120 receives the tuner's output 
receives the captured GPS reference signal s(t)). The down converter 120 
transforms the captured GPS reference signal to an intermediate frequency ("IF") 
reference signal s( t) « To do this, the down converter includes in some embodiments 
20 an IF mixer that converts the frequency of the captured GPS signal to an IF frequency, 
such as 50Mhz. In some embodiments, the down converter also includes one or more 
band pass and amplification stages to filter and amplify the input and/or output of the 
IF mixer. 
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The signal-processing circuit 100 utilizes a down converter so that the A/D 
converter 125 can sample the reference signal at an intermediate frequency as opposed 
to a radio frequency. One of ordinary skill will understand that other embodiments 
can include more than one down converter in their signal-processing circuits. Also, 
5 some embodiments use one or more down converters to convert the GPS reference 
signal to a baseband reference signal, which can then be sampled by the A/D 
converter. 

In some embodiments, the A/D converter's sampling rate is at least twice tiie 
size of the frequency band, while in other embodiments the sampling rate differs from 
10 this rate. The A/D converter 125 samples the IF reference GPS signal s(t ) that it 
receives from the down converter 120, and outputs K-dlgital samples si, sjc of the 
IF GPS signal s(tj. These samples are from only a portion of the received GPS 
signal, and serve as digital GPS reference data that the location determination 
processes use. 

15 These samples are taken at GPS times tj, ...Jx- The corresponding receiver 

clock times are denoted by 7^. . Hence, ^ represent the receiver clock time 
recorded when the first sample of the GPS signal is received at the receiver. The first- 
sample time 7^ represents the approximate time for when the first sample was 
generated because it is measured according to the receiver's clock, which is not 

20 completely synchronized with the GPS time in the embodiments described below. 

The A/D converter also outputs the first-sample time if . In this document, a "tilde" is 
generally placed over a variable (like ^ ) to indicate that the value of the variable is an 

approximation of the actual value of the item represented by the variable. 

13 
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n. ESTIMATING THE LOCATION OF THE GPS RECEIVER 

Figure 2 illustrates a location-detenniiiation process 200 that estimates the 
location of the GPS receiver from the digital samples generated by the signal- 
processing circuit 100. In the embodiments described below, a location-determination 
server separate from the GPS receiver performs this process 200. One of ordinary skill 
will also realize that in other embodiments, the process 200 is performed either 
(1) completely by the GPS receiver, or (2) partially by the GPS receiver and partially 
by another device or computer a location-determination server) in 
communication with the GPS receiver. 

Figure 3 illustrates a location-determination server 300 that performs the 
process 200. In some embodiments, this server is just one computer, while in other 
embodiments several computers form this server. In some embodiments, these several 
computers may be geogr^hically distributed. This server can be a standalone device 
or it can be part of otlier devices. 

As shown in Figure 3, the location-determination server 300 receives digital 
GPS reference data s/, sx from client devices 305 through one or more base stations 
3 1 0. A client device 305 is either a GPS receiver, contains a GPS receiver, or 
.commumcatively connects to a GPS receiver. In addition, each base station detects 
GPS-receiver signal transmissions withia a particular region, and relays this 
information to the location-determination server. One of ordinary skill wiU realize 
that the base station can use a variety of communication architectures and networks to 
relay the signals from the GPS receivers to the location-determination server. In some 
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embodiments, the base station is a cell tower that provides wireless commutiicatioiis 
to the client device 305. 

The location-determination server 300 performs the process 200 for a 
particular GPS receiver whenever it receives the digital GPS reference data Si, sk 
from the particular GPS receiver. In addition to the receiver's reference GPS data, the 
server 300 uses several other data items to perform process 200, in some embodiments 
of the invention. For instance, some embodiments use an approximate location of the 
GPS receiver to compute pseudoranges for several satellites. In the embodiments 
described below, the approximate location is represented in terms of a three- 
dimensional vector 1 , which represents a point in an Earth-centered, Earth-fixed 
coordinate system. Some embodiments specify this approximate location to be within 
±150 Km of the actual GPS receiver location- 
Different embodiments identify the approximate location 7 differently. In the 
embodiments described below, the approximate location is the location of the base 
station 310 that relays the GPS reference data generated by the receiver. Some of 
these embodiments extract the approximate-location information from the signal 
forwarded by the base station. In other embodiments, the server 300 uses the base- 
station's identification to retrieve the base-station location information from a storage 
structure (such as database 315 of Figure 3) that stores the location of all base 
stations. Other embodiments determine the approximate location 7 by measuring 
signal strength between the base station 310 and the client device 305, using time- 
difiference of arrival (TDOA), angle of arrival (AOA), and/or sector of arrival (SOA) 
methods at the base stations. Yet other embodiments identify the approximate 

15 
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location as the last GPS-receivefs location that was recorded within a predetennined 
time interval of the current first-sample time. One of ordinary skill will realize that 
other embodiments use numerous other approaches to identify the approximate OPS- 
receiver location, 

5 In some embodiments, the data-processing server also needs to have access to 

differential GPS data, satellite clock-correction data, navigation bits, and ephemeris 
data. The server can receive such data firom a variety of sources. For instance, as 
shown in Figure 3, the server can receive such data from one or more reference GPS 
receivers 320 through one or more communication networks 325, such as the Internet. 

10 To perfomi the process 200 for a particular GPS receiver, the location- 

determination server 300 initially computes (at 205) the approximate location of all 
GPS satellites at the time that these satellites transmitted the signals resulting in the 
first sample sj of the GPS reference data. This computation is further described below 
in Section II. A by reference to Figure 4. 

15 The process next identifies (at 21 0) N GPS satellites that are over the GPS 

receiver's ^proximate location 1 . This identification is further described below in 
Section n.B by reference to Figures 5 and 6. After identifying the N overhead 
satellites for the approximate location 1 , the process 200 computes (at 215) phase 
offsets for P satellites, where P is a number equal to or less than N, This computation 

20 is further described below in Section II.C by reference to Figure 7. 

For each of the P satellites, the process then computes (at 220) a pseudorange 
firom its computed code phase, the GPS-receiver's approximate location, and the 

16 



wo 02/46786 PCTAJSOl/45828 

satellite's approximate location. As further described below, each satellite's 
approximate location is computed based on the approximate first-sample time ^ , 
which is based on the receiver clock. The computation of the pseudoranges is further 
described below in Section II.D by reference to Figures 8-10. 

The process 200 then uses (at 225) the computed pseudoranges to calculate the 
GPS receiver's estimated location. Some embodiments use triangulation techniques to 
calculate the GPS receiver's estimated location fiom the computed pseudoranges, as 
described below in Section n.E. 

A. Computation of Satellite Locations. 

As mentioned above, the location-determination process 200 initially 
computes flie approximate locations of the satellites at the time that these satellites 
transmitted the signals that resulted in the first sample Sj, The receiver generated this 
first sample at time 7^ . 

Given ephemeris and differential data, a location y' of a satellite i is for all 

practical purposes a deterministic fimction of the GPS time L One set of equations for 
deriving a satellite's location from ephemeris and differential data is provided in Table 
2.3 on page 38 of "Understanding GPS Principles and Applications," by Elliott 
Kaplan, Artech House, 1996. The computed satellite location is typically a three- 
dimensional vector (y € 91^ ) that is defined in an Earth-centered, Earth-fixed 
coordinate system. 
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To compute the location of a satellite j, the process 200 needs to first compute 
the time that this satellite transmitted the signal resulting in the first sample. However, 
the exact transmission time is difiBcuIt to compute without performing signal tracking. 
Hence, the process 200 computes an approximate transmission time and location for 
each satellite, instead of computing the precise values of these parameters. 

More specifically, the exact tratismission time for each satellite i equals the 
actual first-sample time ti minus the actual signal-transit delay r/ for the satellite. Tlie 
signal-transit delay is the time between the satellite's transmission of the signal 
resulting in the first sample sj and the receiver's generation of the first sample. 
Equation (1) illustrates one maimer of computing the exact signal transit delay r/of a 
satellite i. 



In this equation, ^ is the receiver's exact position at first-sample time f/, 3;' (ti-x^ is 
the satellite's precise location at exact transmission time (tj-ri), c is the speed of Ught, 
Tf" is delay due to atmospheric conditions, r''is delay caused by analog processing at 
the receiver. 

The location-determination process 200 can determine atmospheric-delay 
t/' by using differential data. It can also retrieve fcom a storage structure (such as data 

storage 3 15) the processing-delay r'* associated with each GPS receiver. However, it 
is difficult to ascertain (1) the exact first-sample time 0, (2) the receiver's precise 
location £ at exact first-sample time or (3) the satellite's precise location (tj-r^ 




(1) 
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at precise transmission time (f/- 1}). Hence, it is difiBcult to solve the Equation (1) to 
obtain an exact value for the signal-transit delay r/. 

Therefore, in some embodiments, the location-determination process 200 (at 
205) computes for each satellite / the approximate signal transit delay rj of the 
satellite. Equation (2) below illustrates one manner of computing an approximate 
signal-transit delay fj . This equation assumes that (1) the first-sample time is 7^ , 

(2) the receiver is located at the approximate location 1 , and (3) the satellite's 
location at transmission time equals its location at the approximate first-sample time 

r,.H(i,y(i;))=PlM+.,-«' ,2) 

c 

As mentioned above, some embodiments use the base-station tower location as the 
approximate GPS-receiver location 1 , and retrieve the tower location firom a database 
by using the received cell tower identification. 

For each satellite z, this process can thus compute an approximate transmission 
time (7^ - tJ based on the approximate first-sample time 7^ and the satellite's 
approximated signal-transit delay . Accordingly, for each satellite /, the process 
computes the satellite's approximate location (t^-Vf) at the computed approximate 
transmission time (7^ - r J, by using the ephemeris and difTerential data As mentioned 
above, one set of equations for deriving a satellite's location given ephemeris and 
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differential data are provided in Table 2.3 on page 38 of "Understanding GPS 
Principles and Applications," by Elliott Kaplan, Artech House, 1996. 

It can be shown that the approximate satellite location that is computed under 
this approach is no more than few meters off fiom the true satellite location at the time 
5 of transmission. Approximate location errors of less than 150 km will lead to no loss 
in accuracy, while errors in the receive time are less than four meters per millisecond. 
When receiver are part of a network (such as a cellular network), the receiver clocks 
are much more accurate than one nnillisecond. In the remainder of this document, the 
symbol y is used as shorthand for the ap^ximate satellite location y' (t^ - rj. 

10 Figure 4 illustrates a process 400 that uses the above-described approach to 

calculate fee approximate location of each GPS satellite i at the approximate time (7^ - 
fj that the satellite transmitted the signal resulting in the first generated sample. This 
process is used by process 200 at 205, in some embodiments of the invention. 

As shown in Figure 4, the process 400 initially selects (at 405) a first GPS 
15 satellite / from a list of GPS satellites that it has. The process then determines (at 410) 
the signal-processing delay r'' associate with receiver. Some embodiments store in a 
storage structure the processing-^ielay r'' associated with each GPS receiver that can 
be used with the process 200 and the location-detennination server 300. Hence, in 
these embodiments, the process 400 retrieves from the storage structure the 
20 processing-delay t*" associated with the particular GPS receiver at issue. 

Next, for the selected satellite, the process computes (at 415) the atmospheric 

delay r/ . This atmospheric delay can be derived from the received differential GPS 
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data and the approximate receiver location 7 . One manner for deriving atmospheric 
delay at a location from received differential GPS data is disclosed in Chapter 8 of 
"Understanding GPS Principles and Applications," by Elliott Kaplan, Artech House, 
1996. 

5 At 420, the process then computes the selected satellite's approximate location 

y' (T^J aX the approximate first-sample time ij^ . In some embodiments, the process 
computes this location at the time ^ from the ephemeris and differential data by using 
the set of equations provided in Table 2.3 on page 38 of "Understanding GPS 
Principles and Applications," by Elliott Kaplan, Artech House, 1996. 

1 0 After computing the approximate location of the satelUte at time ^ , the 

process 400 computes (at 425) the approximate signal-transit delay t. of the selected 
satellite. This process uses the above-described Equation (2) to compute the 
approximate signal-transit delay from the approximate satellite location y' 

(computed at 420), the signal-processing and atmospheric delays r'' and r,/' 
1 5 (respectively obtained at 41 0 and 4 1 5), and the approximate receiver location I . 

Next, the process computes (at 430) the approximate transmission time (T^-r^) 
for the selected satelUte by subtracting the approximate signal-transit delay 
(computed at 425) from flie approximate first-sample time \ . The process then 
computes (at 440) the selected satellite's approximate location y* at the approximate 
20 transmission time Ot-rJ computed at 43S. Li some embodiments, the process 

computes this location from the ephemeris and differential data by using tbe set of 
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equations provided in Table 2.3 on page 38 of "Understanding GPS Principles and 
Applications," by Elliott Kaplan, Artech House, 1996. 

After computing the selected satellite's location at the approximate 
transmission time, the process determines (at 440) whether the selected satellite is the 
5 last GPS satellite on its GPS-satellite list. If not, the process selects (at 445) another 
GPS satellite from this list, and repeats 415 to 435 in order to determine this newly- 
selected satellite's approximate location. Otherwise, the process end as it has 
computed location of all the GPS satellites at the approximate time they transmitted 
their signals that resulted in the first generated sample. 

10 B. Identification of Overhead Satellites . 

After identifying the approximate location of the GPS satellites, the process 
200 identifies (at 210) the satellites that are currently "overhead." In some 
embodiments, the process 200 identifies the "overhead" satellites by making a 
simplifying assumption that satellites are overhead relative to the client device 305 if 
1 5 and only if they are overhead relative to the approximate location 1 . 

Hence, in these embodiments, the process 200 examines each satellite i and 
determines whether the satellite is over the approximate location 7 . In some 
embodiments, the process specifies a satellite i as an overhead satellite if and only if 
l'l<l'y', (Note that the **prime" denotes vector transposition.) In other words, the 
20 process designates satellite / as an overhead satellite if and only if the inner product of 
the satellite's approximate position vector Si with the approximate location vector 1 is 
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greater than or eqtial to the inner product of the approximate location vector with 
itself. 

This designation criterion essentially determines whether tiie magnitude of the 
projection of a satellite's approximate position vector y onto the approximate 
5 location vector 7 is greater than or equal to the magnitude of the approximate 
location vector £ . If so, the satellite is an overhead satellite. 

Figure 5 pictorially illustrates this designation criterion. This figure presents a 
GPS receiver 505 at the approximate location 7 on Earth 5 10, and two satellites 515 
and 520 that orbit the Earth. The satellite 5 1 5 is not an overhead satellite as it is 

10 beneath the horizon 525 of the location 1 . Accordingly, the magnitude of the 

projection of its position vector ^ onto the approximate location vector t is less 
than the magnitude of the approximate location vector £ . On the other hand, the 
satellite 520 is an overhead satellite as it is above the horizon 525 of the location 7 . 
Hence, the magnitude of the projection of its position vector y^^^ onto the approximate 

15 location vector 7 is greater than the magnitude of the qjproximate location vector I . 

Figure 6 illustrates a process 600 that uses the above-described approach to 
identify the overhead satellites. This process is used by process 200 at 210, in some 
embodiments of the invention. As shown in Figure 6, the process 600 initially 
computes (at 605) the inner product of the approximate location vector I with itself. 

20 The process 600 next selects (at 6 1 0) a GPS satellite from a list of GPS 

satellites that it has. The process then computes (at 615) the inner product of the 
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approximate location vector 7 with the selected satellite's ^proximate location 
vector , which was computed at 205. 

At 620, the process 600 determines whether the inner product generated at 615 
is greater or equal to the inner product generated at 605. If not, the process transitions 
to 630, which will be described below. Otherwise, the process designates (at 625) the 
satellite selected at 610 as an overhead satellite, and then transitions to 630. 

At 630, the process determines whether it has examined all the GPS satellites 
on its list of GPS satellites. If not, the process returns to 610 to select another GPS 
satellite from this list, and repeats the above-described operations in order to 
determine whether the newly-selected satellite is an overhead satellite or not. 

Once the process determines (at 630) that it has examined all the GPS 
satellites, the process ends. Typically, by the time that the process 600 ends, this 
process has identified some integer number of the GPS satellites as overhead 
satellites. 

C. Computation of Phase offsets . 

After identifying the overhead satellites for the approximate location I , the 
process 200 computes (at 215) phase ojffsets for the overhead satellite. One method 
for computing phase ofEsets is described below. However, one of ordinaiy skill will 
realize that a variety of other methods may be used to compute phase of&ets. 

A phase offset for a satellite is the phase difference between (1) the satellite's 
PRN code that is part of the received GPS reference data, and (2) a replica of the 
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satellite's PRN code that is used by process 200. The process 200 computes a 
satellite's phase offset by computing an ambiguity function A^(^, f^) , which is the 
square of the magnitude of the inner product R^ f^) between the received GPS 
reference data with the replica of the satellite's PRN code at various candidate phase 
5 offsets 4> and Doppler-shift values . Equations (3) and (4) provide mathraiatical 
representations of the ambiguity and correlation functions. 

A(lfo)-\Ri(lfo)\'^ (3) 
Ri(lfo) = - ^ - (4) 

In these equations, (1) r^- denotes the GPS time of the h}^ sample, (2) s-j^^ is the A^'' 

1 0 sample after the candidate Doppler-shift value has been applied to it, (3) the candidate 
Doppler-shift vedue/^ represents the sum of the satellite Doppler and the clock 
Doppler values, and (4) x^ft) is a bandpass filtered version of the satellite's C/A code 
replica x^(t) (modulated by the navigation bits) that the process 200 uses. The 
Doppler-shift value is applied to die ^ sample by multiplying this sample by 

1 5 ^-'^'Vr//F*/oXtjr'f,J 

In some embodiments, the search through candidates ((^,^) is exhaustive, 
and the phase of&et of the candidate pair that resulted in the maximimi ambiguity 
value for the satellite is selected as the phase offset of the satellite. As described 
below by reference to Figure 7, some embodiments take this selection a step fiirthCT. 
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These embodiments use an ambiguity tibreshold ratio to detect useful phase oflfeets, 
and disregard satellites without useful phase offsets. For instance, some embodiments 
reqioire a satellite's maximimi ambiguity to exceed the average ambiguity by a 
constant ratio in order for its phase ofifset to be considered valid. In other 
embodiments, a variety of threshold tests can be used, such as ratio between highest 
ambiguity peak and next highest ambiguity peak. 

Also, rather than computing the inner-product function Rf($,fQ) as a sum of 
a series of multiplications for each ^ and , some embodiments view as a 

correlation. In one embodiment, the correlation theorem and fast Fourier transforms 
("FFT") are used to improve computation time. 

Figure 7 presents a process 700 for computing useful phase offsets for several 
of the overhead satellites. The process 200 performs (at 215) the process 700 after it 
has identified the overhead satelUtes. The process 700 initially (at 705) selects an 
overhead satellite and defines this satellite's maximum ambiguity value to be equal 
toO. 

Next, the process 700 enumerates (at 710) a number of candidate phase ofi&ets 

^ *s and Doppler-shifl values 's for the selected satellite. The process then selects 

(at 715) a candidate phase of&et ^ and Doppler-shift value , and computes (at 

720) the ambiguity value for the selected candidate pair by using Equations (3) and (4) 
above. 

The process next determines (at 725) whether the ambiguity value computed at 
720 is greater than the maximum ambiguity value recorded for the selected satellite 
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thus far. If not, the process transitions to 735, which will be described below. If so, 
the process (at 730) sets the maximum ambiguity value equal to the value computed at 
720, and stores the satellite's phase offset as the best phase offset thus fer. The 
process then transitions to 735. 

At 735, the process determines whether it has examined all the candidate pairs. 
If not, the process returns to 715 to select another candidate pair. Otherwise, the 
process determines (at 740) whether the selected satellite's maximum ambiguity value 
exceeds a threshold ambiguity value. If so, the process adds (at 745) the selected 
satellite to the list of satellite's for which a pseudorange needs to be computed, and 
records the best phase offset stored at 730 for this satellite as its phase offset. 

From 745, the process transitions to 750, The process also transitions to 750 if 
it detemiines (at 740) that the selected satellite's maximum ambiguity value does not 
exceed the threshold ambiguity value. At 745, the process determines whether it has 
examined all the overhead satellites. If not, the process returns to 705 to select another 
satellite. Otherwise, the process ends. 

D. Comptttation of Pseudoranges. 

After computing phase offsets for several overhead satellites, the process 200 
then uses (at 220) the receiver's approximate location and the computed phase of&ets 
to compute pseudoranges for the set of satellites. To compute a satellite's 
pseudorange, the process 200 calculates the length of all or most of the PRN codes 
transmitted by the satellite from the time that it transmitted the signal resulting in the 
first sample si to the time the receiver generated the first sample sj. 
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As mentioned above, signal-transit delay refers to this tintie duration (i.e., it is 
the time between the satellite's transmission of the signal resulting in the &st sample 
sj and the receiver's g^eration of the first sample si). During the satellite's signal- 
transit delay, satellite i transmits several PRN codes. 

S Figure 8 Ulustrates that these PRN codes can be grouped in three categories. 

The first category is the PRN code segment 805 that contains the signal resulting in 
the first sample sj. In the discussion below, this PRN code segment is referred to as 
the first PRN code segment. The first PRN code segment 805 is said to be a complete 
or practically complete segment when the signal resulting in the first sample is part of 
10 the first PRN-code chip, and it is said to be a partial segment when the signal resulting 
in the first sample is part of a chip after the first PRN-code chip. 

The second category includes the full PRN codes 810 tliat the satelUte 
transmitted after the first PRN code but before the first-sample time. In the discussion 
below, these PRN codes are referred to as the full PRN codes. The third categoiy is 

15 the PRN code 8 1 5 that the satellite is transmitting when the receiver generates the first 
sample. In the discussion below, this PRN code segment is referred to as the last PRN 
code segment 815. The last PRN code segment 815 is said to be a complete or 
practically complete segment when at the first-sample time the sateUite is transmitting 
the last chip of this code, and it is said to be a partial segment when at the first-sample 

20 time the satellite is transmitting a chip before the last chip of this code 

Equation (5) below illustrates one manner of computing the estimated distance 
between the GPS receiver and the satellite by summing the length first, last, and 
full PRN code segments after accounting for differential corrections and the satellite- 
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clock's error. In other words. Equation (5) provides the mathematical representation 
of a pseudorange p^ for the satellite. 



In this equation, (1) is the length of the first code segment, (2) is the number of 
ftdl PRN codes, (3) / is the full-PRN-code lengfli, which is about 300 Km, (4) is 
the sateUite's transmission phase, (5) c is the speed of light, (6) cv, is the length of 
the last PRN code without adjustment for the satellite clock error, (7) is the 
satellite-clock-error correction factor, and (8) is tfie atmospheric-error-correction 
factor. 

Accordingly, in Equation (5), (1) the first term ( /?, ) accounts for the length of 
the first code segment, (2) the second term (m^/) accoimts for the length of the full 
PRN codes, (3) the third term ( cv. + rf/ ) accounts for the length of the last PRN code 
after adjustment for satellite-clock error, and (4) the fourth term (d, ) accounts for 
atmospheric delay. 

Figure 9 illustrates a process 900 for computing a satellite's pseudorange 
according to Equation (6) below. 



Equation (6) provides the mathematical representation for computing another 
pseudorange for the satellite. Equation (6) is similar to Equation (5) except tiiat 



Pt^-Pf-^mJ + fcvt-^d'j + dt 



(5) 



Pf = Pi'{-mfl'hd/' -hdf 



(6) 
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it does not include cv, , which is the Iast-PRN-code*s length wifliout adjustment for 
satellite clock error. As further described below, the process 900 does not compirte 
cVf , but rather leaves it to a subsequent triangulation operation to account for this 
tenn for all the P satellites together. 



The length Pf of the first code segment 805 of satellite i is the code length 
associated with the phase offset of this satellite. Accordingly, this code length is 
referred to as phase length in the discussion below. Figure 10 illustrates the 
relationship between the phase length and phase offset of satellite /. This figure 
presents three codes. The first code 1005 is the ^ satellite's PRN code in the received 
signal. (In this figure, the illustration of the received PKN code disregards 
degradation and distortion in the received signal.) 

The second code is a replica 1010 of the satellite's PRN code that the 
location-determinadon process uses. Figure 10 illustrates the length of complete code 
replica. The third code is a version 1015 of the replica that has been shifted by a 
phase offset to perfectly match the received code 1005. As shown in Figure 10, the 
phase length associated with the phase offset equals the length of the first code 
segment in the received signal. 

When the phase offset is represented as a nonnalized-time value between 0 
to 1, the /'^ satellite's phase length is represented by Equation (7) below: 



1- 



Computing the length of the first PRN code segment 



Pi = <fiiC . 



(7) 
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The phase lengfli can be expressed differently. For instance, it can be expressed 
differently when the phase ofiEset is represented differently. 



To compute the length of the full PRN codes 810 that the f satellite 
transmitted after the first PRN code 805 but before the first-sample time, some 
embodiments first compute the number of these complete codes. Equation (8) 
provides a mathematical representation for calculating this number. 



As illustrated by this equation, the number of full PRN codes can be estimated based 
on the i^^ satellite's approximate location y' , the GPS-receiver*s approximate location 
i , the satellite's computed phase length , and the PRN-code-length /. 

Specifically, this number can be estimated by (1) subtracting the satellite's 
approximate location firom the GPS-receiver*s approximate location to obtain an 
approximate distance between the satellite and receiver, (2) subtracting the phase 
length plus distance correction factor d' from this distance, (3) dividing the resulting 
value by the PRN-code-length /, and (4) rounding the result of the division to the 
closest integer value. The calculation of the distance correction factor is 
discussed in the following subsection. Once the number of full PRN codes is 
estimated, the length of these codes can be obtained by multiplying this number by the 
PRN-code-length /, as mentioned above. 



2. 



Computing the length of full PRN code segments 
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3. 



Computing the length of tibe last PRN code segment and 
adjusting for satellite-clock error. 



So long as the GPS-satellite clocks are all synchronized by accoxmting for the 
satellite-clock errors, the length of the last code segment 8 1 S is equal for all P 
satellites. Hence, as mentioned above, the process 900 below sets the unadjusted last- 
PRN-code length (cv^ ) to zero, and leaves it to a subsequent triangulation operation to 
account for this term for all the P satellites together. 

As mentioned above. Equation (6) includes a distance correction factor d^ to 

account for the satellite-clock errors. Equation (9) provides one formula for 
calculating this error correction factor. 



In this equation, is the satellite-clock error, which in some embodiments is 

retrieved by the location-determination process fix>m a GPS receiver 320 through the 
communication networks 325. 



As mentioned above. Equation (6) includes a differential correction factor . 
This correction factor is to account for signal-path delays due to atmoi^heric 
conditions and other effects. This diJSerential correction factor can be obtained from a 
variety of sources, including third-party vendors, publicly available data via the 
communication network 325 (e.g, www.ngs.noaa.gov/ CORS/cors-data.html). readings 
of an indqsendent GPS receiver and other sources. 




(9) 



4. 



Adjusting for differential corrections. 
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5. Pseudorange-computing process 900, 

The pseudorange-computing process 900 of Figure 9 will now be described. 
This process initially selects (at 905) one of the P satellites that the process 700 
selected for computing their pseudoranges. Next, the process uses (at 910) the above- 
described Equation (7) to compute the selected satellite's phase length from its 
phase offset ^^ , which was computed by process 700. 

The process then xises (at 915) Equation (9) to compute distance correction 
factor rf/ due to the satellite-clock error. As mentioned above, some embodiments 
retrieve the satellite-clock error from a GPS receiver 320 through the 

communication networks 325. At 920, the process uses Equation (8) to compute the 
number of the full PRN codes 810 that the selected satellite transmitted after the first 
PRN code 805 but before the first-sample time. The process 900 next computes 
(at 925) the length of these full PRN codes by multiplying the number computed at 
91 5 by the PRN-code-length /, which is about 300 Km. At 930, the process then 
identifies an atmospheric-delay correction factor . 

The process then computes (at 935) the selected satellite's pseudorange by • 
summing the values computed at 910, 920, 925, and 930. The process 900 then 
determines (at 940) whether it has generated pseudoranges for all the P satellites 
identified by process 700. If not, the process returns to 905 to select the next satellite, 
and repeat 910-935 to compute the newly-selected satellite's pseudorange. Otherwise, 
the process ends. 
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Triangulation of Pseudoranges to Determine LocatioiL. 

After computing pseudoranges, the process 200 then uses (at 225) the 
computed pseudoranges to calculate the GPS receiver's estimated location. Some 
embodiments use triangulation techniques to calculate the GPS receiver's estimated 
location from the computed pseudoranges. 

There are multiple triangulation approaches, such as those discussed on pages 
44-48 of "Understanding GPS Principles and Applications,*' by Elliot Kaplan, Artech 
House, 1996. One skilled in the art will realize that any one of these approaches can 
readily be used in conjunction with the invention to compute receiver location from 
the computed pseudoranges. 

In some embodiments, the triangulation algorithm uses four satellites to 
compensate for a common uncertainty in the satellite distances, instead of 
compensating for the receiver clock error, which is the traditional GPS interpretation 
for using four satellites during triangulations. As discussed in Section n,D and Figure 
8, tills uncertainly is the last PRN Codes 815. TTiis term is common to all the 
satellites, and thus, by including four satellites, the triangulation algorithm can solve 
for this term. Hie receiver clock error and the uncertainty of &e last PRN Code are 
equivalent manifestations of the same entity. 

One of ordinary skill will realize fbat the above-described embodiments have 
several advantages. For instance, these embodiments do not perform signal tracking 
to compute the pseudoranges of satellites. Accordingly, these embodiments are faster 



34 



wo 02/46786 PCT/USO 1/45828 

than some of the prior techniques. In addition, these embodiments do not need several 
seconds of GPS data, as they do not extract the satellite time stamps. 

While the invention has been described with reference to numerous specific 
details, one of ordinary skill in the art will recognize that the invention can be 
embodied in other specific forms without departing fi-om the spirit of the invention. 
For instance, even though the embodiments described above use the commercial GPS 
frequency, other embodiments use the military GPS firequency, while others use both 
the military and commercial fi-equencies. Thus, one of ordinary skill in the art will 
tmderstand that the invention is not to be limited by the foregoing illustrative details, 
but rather is to be defined by the appended claims. 
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CLAIMS 

We claim: 

1 . A global positioning method comprising: 

a) using a global positioning system ("GPS") receiver to receive a 

GPS signal; 

b) computing a pseudorange for a GPS satellite without 
performing signal tracking on the received GPS signal; and 

c) using the computed pseudorange to identify an estunated 
location of the GPS receiver. 

2. The global positioning method of claim 1 further comprising: 

a) computing pseudoranges for at least two other GPS satellites 
without performing signal tracking on the received GPS signal; and 

b) using the computed pseudoranges to perform triangulation to 
identify the estimated location of the GPS receiver. 

3 . A method of computing a pseudorange for a GPS receiver and a GPS 
satellite, vsiierein the GPS receiver receives a GPS signal and the GPS satellite 
repeatedly transmits a code, the method comprising 

a) identifying a phase offset between the received GPS signal and 
a replica of the satellite's code, wherein the phase oflfeet corresponds to a particular 
segment of the code; 
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b) identifying an approximate GPS-receiver location; 

c) identifying an approximate time for the received GPS signal; 

d) computing the pseudorange for the GPS satellite based on the 
computed phase ofiset, the approximate GPS-receiver location, and the approximate 
time. 

4. The method of claim 3, wherein the approximate GPS-receiver location 
is within 1 50 Km of the actual receiver location. 

5. The method of claim 3, wherein the receiver has a clock that maintains 
time, the method further comprising maintaining the receiver-clock time to be within 
a predetermined time interval of GPS time. 

6. The method of claim 5, wherein medntaining the receiver-clock time 
comprises synchronizing the receiver-clock time to match time generated by a 
reference clock. 

7. The method of claim 6, wherein synchronizing the receiver-clock time 
comprises periodically synchronizing the receiver-clock time to match time generated 
by the reference clock. 

8. The method of claim 5, wherein the predetermined time interval is 1 
millisecond. 

9. A method of computing a pseudorange for a receiver and a transmitter, 
wherein the receiver receives a reference signal starting at a particular time and the 
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transmittef trazismits. a signal with a periodically-repeating component, the meOiod 
comprising 

a) identifying a phase offset between the received reference signal 
and a replica of the transmitter signal's periodically-repeating component, wherem the 
phase of^t corresponds to a particular segment of the periodically-repeating 
component; 

b) computmg the length of the particular segment; 

c) computing the length of complete periodically-repeating 
.componen,ts that the particular satellite has transmitted by the particular time after 
transmitting the particular segment; 

d) calculating the pseudorange based on the computed lengths. 

10. The method of clahn 9, wherein calculating the pseudoranges 
comprises summing tlie computed lengths. 

1 1 . The method of claun 9, wherein identifying the phase offset includes 
correlating the received signal with a replica of the transmitter signal's periodically- 
repeating component at various phase offsets. 

1 2; The method of claim 9, wherein computing the length of complete 
periodically-repeating components comprises: 

a) calculating the number of complete periodically-repeating 
components that the particular satellite transmitted after the particular segment by the 
particular time; and 
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b) multiplying the calculated number by the length of the complete 
periodically-repeating component. 

1 3 . The method of claim 9, wherein the transmitter is part of a GPS 
satellite and the receiver is a GPS receiver, 

14. The method of claim 1 3 , wherein the satellite includes a clock and the 
clock differs from GPS clock by some error amount, the method further comprising: 

a) identifying the satellite's clock error amount; 

b) genemting an error-correction factor; 

c) taking the error-correction factor into account in calculating the 

pseudorange. 

1 5 . The method of claim 1 4, 

a) wherein generating an error-correction factor comprises 
multiplying the error amount by the speed of light; 

b) wherein taking the error-correction factor into account 
comprises adding the error-correction &ctor to the calculated pseudorange. 

' 1 6. The metihod of claim 1 3, wherein calculating the pseudorange further 
comprises: 

a) calculating an atmospheric-delay-correction &ctor for the 

satellite; 
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b) taking the correction factor into account in calculating the 

pseudorange. 

17. The method of claim 16, wherein taking the correction factor into 
accoimt comprises adding the correction factor to the calculated pseudorange. 

5 18. For a global positioning system ("GPS") that includes a plurality of 

GPS satellites and a GPS-receiver, a method of identifying an estimated location of a 
GPS receiver, wherein each particular GPS satellite transmits a signal that includes a 
periodically-repeating code unique to the particular GPS satellite, wherein the GPS 
receiver receives a signal starting at a particular time, the method comprising 

10 a) for each particular GPS satellite in a set of the GPS satellites, 

1 ) determining a phase offset between the received signal 
and a replica of the satellite's code, wherein the phase offset corresponds to a 
particular segment of the particular satellite's code; 

2) computing the length of the code segment; 

15 3) computing the length of complete codes that the 

particular satellite transmitted after the code segment by the particular time; and 

4) calculating a pseudorange for the particular satellite 
based on the computed lengths; and 

b) identifying the estimated location of the receiver by using the 
20 computed pseudoranges. 
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19. Tlie method of claim 18, wherein calculating the pseudorange for each 
particular satellite comprises summing the computed lengths for the particular 
satellite. 

20. The method of claim 1 8, wherein identifying the phase offeet for each 
particular satellite includes correlating the received GPS signal with a replica of the 
particular satellite's code at various phase of&ets. 

21 . The method of claim 1 8, wherein computing the length of complete 
codes for each particular satellite comprises: 

a) calculating the nxmiber of complete codes that the particular 
satellite transmitted after the code segment by the particular time; and 

b) multiplying the calculated number of complete codes by the 
length of the satellite's code. 

22. The method of claim 18, wherein each satellite includes a clock and the 
clock differs fix)m GPS clock by some error amount, the method further comprising: 

for each particular satellite, 

a) identifying the particular satellite'^ clock error amount; 

b) generating an error-coirection factor; 

c) taking the error-correction factor into accoimt in calculating the 
pseudorange for the particular satellite. 

23. The method of claim 22, 
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a) 



wherein generating an error-correction &ctor comprises 



multiplying the error amount by the speed of light; 

b) wherein taking the error-correction factor into account 
comprises adding the error-correction factor to the calculated pseudorange for the 
particular satellite. 

24. The method of claim 1 8, wherein calculating the pseudorange for each 
particular satellite further comprises: 



particular satellite; 

b) taking the correction factor into account in calculating the 
pseudorange for the particular satellite. 

25. The method of claim 24, wherein taking the correction factor into 
account comprises adding the correction factor to the calculated pseudorange for the 
particular satellite. 



a) 



calculating an atmospheiic-delay-conection factor for the 
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